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O ' Abstract 

<N , 

, An open question remains whether Ultraluminous X-ray Sources (ULXs) really contain intermediate- 

Q ■ mass black holes (IMBHs) . We carefully investigated the XMM- Newton EPIC spectra of the four ULXs that 

were claimed to be strong candidates of IMBHs by several authors. We first tried fitting by the standard 
spectral model of disk blackbody (DBB) + power-law (PL), finding good fits to all of the data , in agreement 
■"sj" ' with others. We, however, found that the PL component dominates the DBB component at ~ 0.3 to 10 

04 , keV. Thus, the black hole parameters derived solely from the minor DBB component are questionable. 

Next, we tried to fit the same data by the "p-free disk model" without the PL component, assuming the 
effective temperature profile of T^b oc where r is the disk radius. Interestingly, in spite of one less 
' free model parameters, we obtained similarly good fits with much higher innermost disk temperatures, 

1.8 < kTin < 3.2 keV. More importantly, we obtained p ^ 0.5, just the value predicted by the slim (super- 
. critical) disk theory, rather than p = 0.75 that is expected from the standard disk model. The estimated 

0^ ' black hole masses from the p-free disk model are much smaller; M ^ 40Mq. Furthermore, we applied a 

, more sophisticated slim disk model by Kawaguchi (2003, ApJ, 593,69), and obtained good fits with roughly 

consistent black hole masses. We thus conclude that the central engines of these ULXs are super-critical 
accretion flows to stellar-mass black holes. 
• Key words: accretion, accretion disks — black hole physics — X-rays; individuals (NGC 5204 X-1, 

P-t\ NGC 4559 X-7, NGC 4559 X-IO, NGC 1313 X-2) — X-rays: stars 

o 

' 1. Introduction are single compact objects powered by accretion flow. 

, Then, how can we understand the extremely large lumi- 

X-ray data analysis has become one of the most widely nosities exceeding the Eddington luminosity for a mass of 

studied subjects in these days, and has made great con- ~10 A'/q? 

tributions in the progress of astrophysics. This is particu- A possibility that ULXs are only bright toward our line 

?H I larly true in the field of astrophysical black holes. Thanks of sights (beaming model) was proposed (e.g. King et al. 

to rapid improvement of the sensitivity of X-ray detectors, 2001, Kording et al. 2002, Reynolds et al. 1997), but is 

new types of X-ray sources have been discovered. Some of not generally accepted now. In fact, Wang (2002) has 

them seem to be black holes which were not known before, shown the existence of a bright nebula surrounding the 

Recent X-ray observations on many nearby spiral galax- ULX M81 X-9, and some ULXs also exhibit extended op- 

ies showed that the X-ray radiation consists of emission tical nebulae (PakuU, Mirioni 2002). There is an evidence 

from several discrete X-ray sources, such as hot gaseous that some ULXs are associated with a spreading wave of 

media, active galactic nucleus (AGN) and the point-like star formation, which supports the idea that many ULXs, 

off-center X-ray sources whose X-ray luminosities signif- in particular those located in star forming galaxies, are 

icantly exceed the Eddington luminosity of a neutron high mass X-ray binaries (King 2004). These facts indi- 

star (Fabbiano 1989). The bright point-like off-center cate that the ULX surrounding nebulae are powered by 

X-ray sources which cannot be identified with young the central high mass X-ray sources, and that strong X- 

supernova remnants are known as ultraluminous X-ray ray beaming toward our direction is unlikely, 

sources (ULXs; e.g., Makishima et al. 2000; Ptak, Colbert Then, the question is, how massive the central com- 

2004). The ULXs, despite the fact of being very lumi- pact objects are. At present, there are basically two 

nous (ix ^ 10^^~^^ergs~^). cannot be explained as col- distinct lines of thoughts; sub-Eddington accretion onto 

lections of many sources with each luminosity less than an intermediate (several hundreds Mq) mass black hole 

the Eddington limit, since many of the ULXs exhibit sig- (IMBH), and super- Eddington (super-critical) accretion 

nificant time variabilities (e.g., Makishima et al. 2000). onto a stellar-mass black hole. The success in fitting 

Therefore, one reasonable assumption is that the ULXs several ULXs spectra with a multicolor disk blackbody 
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(DBB) and a power-law (PL) model supports the idea of 
IMBHs for the ULXs (Miller et al. 2003, Miller et al. 2004; 
Cropper et al. 2004, Roberts et al. 2005), since the cool 
inner disk temperatures and the large innermost radii, ob- 
tained through model fitting, suggests a black hole mass 
within the IMBH range. The other line of thoughts, 
supporting a super-Eddington accretion flow, is the slim 
disk model with stellar mass black holes (Watarai et al. 
2001; Kawaguchi 2003, Ebisawa et al. 2003; Okajima et 
al. 2006). The present article concerns with the super- 
Eddington thick accretion disks introduced in the early 
1980s by the Warsaw group and their collaborators (e.g. 
Abramowicz et al. 1978; Abramowicz et al. 1980) and with 
the slim accretion disks, introduced in the late 1980s by 
the Warsaw and the Kyoto groups (Abramowicz et al. 
1988, Abramowicz et al. 1989). In these papers, the very 
possibility of the super-Eddington beaming in the thick 
accretion disk funnels, crucial for arguments presented 
here, was explicitly recognized, discussed and stressed. 

The IMBH notion faces a serious problem that 
the formation of IMBH itself remains unsettled (e.g., 
Madhusudhan et al. 2006). In contrast, the alternative 
possibility of the super-critical accretion has been poorly 
investigated so far due mainly to our limited knowledge 
about properties of the super-critical flow. The situation 
has been remarkably improved in the past few years, how- 
ever, since the basic tools for investigating high-luminosity 
accreting systems are now available. In the present pa- 
per, we investigate the XMM-Newton EPIC spectra of 
the four ULXs which had been claimed as strong IMBH 
candidates, using theoretical super-critical models to see 
whether the super-critical accretion takes place in these 
ULXs or not. Our final goal is to determine the mass of 
the central black holes and settle down the issue regarding 
the origin of ULXs. 

We applied accretion disk spectral models to the ob- 
served ULX energy spectra, constrain the innermost disk 
radius, and estimate the black hole mass. A very similar 
method as the one used in the present paper was success- 
fully applied to Galactic black hole binaries (e.g., Ebisawa 
et al. 1991; Ebisawa et al. 1993; Gierlihski et al. 2001). 
Several authors have recently done elaborate mass and 
angular momentum estimates by spectral fitting for mi- 
croquasars, i.e. by a very similar method as the one used 
in the present paper (Shafee et al. 2006; McClintock et al. 
2006; Davis et al. 2006). In the present study, however, 
we focus on the mass estimation and do not attempt to 
estimate the spin of the black hole. 

The mass of the ULXs can also be accurately estimated 
by a different method than that described in our arti- 
cle. For this other independent estimate, double peak 
QPOs in the 3:2 resonance should be detected in ULXs 
(Abramowicz et al. 2004). This is because these QPOs 
scale inversely with the mass of the source, which was 
first shown for microquasars by McClintock and Remillard 
(2003), and more recently for low-mass Seyfert galaxies 
by Lachowicz et al. (2006). There is also an indication 
that the same scaling is true for a (possible) detection 
of double peak 3:2 QPOs in SgrA* (Torok 2005), while 



for the ULXs, in general, it has not been detected yet. 
According to Mucciarelli et al. (2006), the only ULX 
where a QPO has been discovered at present is M82 X- 

1. Scaling the frequency inversely to the BH mass, the 
observed QPO frequency range (from a few tens to a few 
hundreds mHz) would yield a black hole mass anywhere in 
the interval from a few tens to a few thousands solar mass 
(Strohmayer, Mushotzsky 2003; MucciareUi et al. 2006). 

The Plan of this paper is as follows: We first describe 
the X-ray data which we used in the present study and our 
methods of fitting in Section 2. In Section 3, we then give 
the results of fitting, together with their implications on 
the black hole mass and mass accretion rates. Section 4 
is devoted to discussion, and the final Section 5 concludes 
the paper. 

2. Data and Fitting Methods 

2.1. Data 

We perform spectral analysis of four ULXs by using 
XMM-Newton observational data. The data were ex- 
tracted using the XMM-Newton SAS version 6.5.0 tools. 
Detail extraction process will be explained separately for 
each object. Response files were made using the SAS tools 
rmfgen and arfgen for all the data. 

2.1.1. NGC 5204 X-1 

NGC 5204 X-1 is located ~ 0.3 kpc from the center of a 
nearby, D = 4.8 Mpc, Magellanic-type galaxy (Roberts, 
Warwick 2000). The typical X-ray luminosity of this 
source is of the order of 2 6 x lO'^^erg s~^ (0.5 - 8 keV) 
(Roberts et al. 2005). 

NGC 5204 X-1 data for our analysis were obtained on 
2003 January 6 (observation ID 0142770101). Following 
Roberts et al. (2005), we set flag = and pattern < 4 
for pn data , while #XMMEA_EM flag and pattern < 12 
were used for the MOS data. No time filters were applied 
to the data since the background is less than 10 count 
s~^ in the pn detector. Source spectra and light curves 
were extracted in 36 — arcsec radius circle centered on the 
source and 44 — arcsec radius region nearby was used to 
produce background data. 

2.1.2. NGC 4559 X-7 and X-10 

NGC 4559 X-7 and X-10 are two brightest sources in 
NGC 4559 (D = 9.69 Mpc). X-7 was located at the outer 
spiral arms, while X-10 was about 0.3 kpc from the op- 
tical nucleus of NGC 4559. They are also known as XI 
and X4 respectively in Roberts and Warwick (2000), and 
1X065 and 1X066 in Colbert and Ptak (2002) respectively 
(Cropper et al. 2004). 

The data were obtained on 2003 May 27 (observation 
ID 0152170501) for both X-7 and X-10. However, X-10 
was only seen in EPIC-pn since the EPIC-MOS cameras 
were operated in small window mode. Following Cropper 
et al. (2004), we used 30 — arcsec radius centered to the 
source for X-7, while for X-10 we used 20 — arcsec ra- 
dius with an exclusion of 6 — arcsec radius centered on 
the nucleus of the galaxy (12''35'"57."64; +27°57'35".8, 
J2000.0). Source free region on the same CCD was se- 
lected as background for the EPIC-pn data while region 
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as close as possible to the source were selected for EPIC- 
MOS X-7 data. The MOS data were filtered using pattern 
< 12 and #XMMEA_EM flag, while the pn data were fil- 
tered using pattern < 4 and #XMMEA_EP flag. 
2.1.3. NGC 1313 X-2 

NGC f313 X-2 is one example of ULX sources in a 
nearby spiral galaxy, NGC 13f3 [D = 3.7 Mpc). This 
ULX is located at approximately 8 kpc from the photo- 
metric center of the galaxy (Colbert et al. 1995). The 
data for our analysis were obtained on 2000 October 17 
(observation ID 0106860101). 

We performed the standard procedures but somehow we 
failed to create response files for the MOS data, leaving 
us only the pn data to be analyzed. We set flag = and 
pattern < 4 for the pn-data. Following Miller et al. (2003), 
we used 24-arcsec radius centered to the source to obtain 
source spectra, while background counts were extracted 
in an annulus between 24-arcsec and 30-arcsec. 

2.2. Fitting Methods 

The source spectra were grouped to a minimum of 20 
count/bin before fitting using XSPEC version 11.3. We 
first try the standard spectral model of disk blackbody 
(DBB) -I- power-law (PL). We use the results as compar- 
ison with other papers. Next, we tried to fit the data 
by the "p-frcc disk model", assuming the effective tem- 
perature profile of Toff oc r"^, where r is the disk radius 
(Mineshige et al. 1994). The p-free model is a potentially 
useful spectral model, since despite its simpleness it is very 
powerful for discriminating slim disk with p ^ 0.5 from the 
standard disk with p = 0.75 (Watarai ct al. 2000). 

We also tried to fit the data with a more sophisticated 
slim disk model calculated by Kawaguchi (2003), who cal- 
culated more realistic emergent spectra of super-critical 
flow based on the slim disk model, taking into account rel- 
ativistic and Compton-scattering effects. We will show the 
results separately apart from those two mentioned above, 
to avoid muddle, and we will give a brief comment on 
these results. 

The effect of Galactic absorption was taken into ac- 
count by using wabs model (in XSPEC), where we fixed 
the value at 1.5 x lO^^cm-^ for NGC 5204 X-1 and NGC 
4559 X-7 and X-10, and 3.9 x lO^^cm'^ for NGC 1313 X-2 
from Dickey and Lockman (1990). The effect of absorp- 
tion external to our galaxy was also considered and again 
we used wabs model but we let this value free throughout 
the fitting. 

3. Spectral Analysis 

The results of spectral analysis for each ULX investi- 
gated in this paper will be presented separately for three 
distinct spectral models: DBB-t-PL models (in sections 3.1 
and 3.2), p-free models (in section 3.3), and Kawaguchi's 
models (in section 3.4). In addition, we summarize the 
fitting results in Table 1. 



3.1. DBB+PL model 

Wc first present our results of fitting based on the con- 
ventional approach; i.e., by using the DBB-I-PL model. 
As will be explicitly demonstrated below, we will basi- 
cally confirm the previous results supporting the IMBH 
hypothesis. 

3.1.1. NGC 5204 X-1 

We fit the MOS data over 0.3 - 8 keV and the pn data 
over 0.3 - 10 keV. We obtained a good fit to the data with 
DBB+PL model (see Fig.l.). The reduced chi-squared is 
1.05 with the inner disk temperature, fcTin ~ 0.25 ±0.03 
keV, and photon index, F = 1.92 ± 0.06. The low temper- 
ature obtained from the fitting apparently supports the 
IMBH interpretation of ULXs, in agreement with Roberts 
et al. (2005). The flux obtained from the fitting by using 
this model is, (0.3 - 10 keV)= 1.66 x lO-^^^rg cm-^s"!. 
These are in reasonable agreement with the previous re- 
sults by Roberts et al. (2005); that is, fcTin = 0.21 ±0.03 
keV, and F = 1.97 ±0.07. 

3.1.2. NGC 4559 X-7 and X-10 

We next fit both the MOS and the pn data for NGC 
4559 X-7 and the pn data for NGC 4559 X-10 over 0.3 - 10 
keV. As before, we fit the data with the DBB+PL model 
and we obtained reduced chi-squared of 0.96 and 0.92 for 
X-7 and X-10, respectively. The inner disk temperature, 
fcTin = 0.17 + 0.09 keV, and F = 2.17 ±0.05, were obtained 
for NGC 4559 X-7, while fcTin = 0.48 ±0.25 keV and F = 
1.94 ±0.13 were obtained for NGC 4559 X-10. However, 
for the case of X-10, we could actually obtained a good 
fit by the PL model alone, giving reduced chi-squared of 
0.92 and F== 1.99 ±0.03. 

For reference, the previous fits by Cropper (2004) ob- 
tained kTin = 0.148 ± 0.006 keV, and F = 2.23 ± 0.05 for 
X-7, in good agreement with our results. As for X-10, 
Cropper et al. (2004) only shows the fitting result with 
single component power-law. 

3.1.3. NGC1313 X-2 

We fit the pn data of NGC 1313 X-2 over 0.2 - 10 keV 
with the DBB+PL model and found a good fit with re- 
duced chi-squared of 0.96. As in previous cases, we obtain 
low inner disk temperature, fcTjn = 0.27 ± 0.04 keV, and 
F = 2.01 ±0.13. 

For comparison. Miller et al. (2003) found fcTin = 0.16± 
0.16 keV, and F = 2.3±0.2, for NGC 1313 X-2, from fitting 
the MOS data with the DBB+PL model. The difference of 
the detectors is likely the origin of the small discrepancies. 

3.2. Interpretation based on DBB+PL model 

3.2.1. Basic methodology to derive black-hole mass 

When we obtain a good fit to the observed spectra by 
the DBB model, we can easily estimate the black-hole 
mass and the Eddington ratio, L/L-e (with L and Le 
being disk luminosity and the Eddington luminosity, re- 
spectively), based on the standard disk theory. The basic 
methodology is summarized in Makishima et al. (2000). 

Following Makishima ct al.(2000), the bolometric lumi- 
nosity of an optically thick accretion disk can be written 
as 




Fig. 1. Best fitting spectra and the fit residuals with the DBB+PL model (left panels) and p-free disk model (right panels) for 
NGC 5204 X-1, NGC1313 X-2, NGC 4559 X-7 and X-10. Open circles, open squares and dots represent mosl, mos2, and pn data, 
respectively. Dashed lines in the left panels represent spectral components. 



ibol = 2^i^Vboi(cosi)-i (1) 

with i is the inchnation of the disk (i = corresponds 
to face-on geometry) and D is the distance. This Lboi is 
related to the maximum disk color temperature Tin and 
the innermost disk radius Rin as 

ibol = MRin/0^<y{Tin/^f (2) 

where (Shimura, Takahara 1995) is the ratio of the 

color temperature to the effective temperature, or spectral 
hardening factor, and ^ — 0.412 is correction factor reflect- 
ing the fact that Tin occurs at somewhat larger than Rm- 
Hence, the innermost disk radius is 

We may identify with the radius of the last stable 
Keplerian orbit. Thus, we may in general write 

M 

Rin = 3/3i?s = 8.86/3(— )km (4) 

Mq 



where Rs is the Schwarzschild radius (= 2GM/c'^), by 
which we can determine the black hole mass. Note that 
1/6 < P < 1, /3=l/6 for extremely rotating Kerr black 
hole and /? = 1 for Schwarzschild black hole. 

Here, we use 0.3 - 10 keV flux, /[o.3-iokeV]i from which 
we derive 0.3 - 10 keV luminosity, i[o.3-iokeV] • We assume 
Lhoi ~ -^[0.3-iokoV] because all significant energy range 
from 0.3 to 10 keV. Therefore, using Tin and /[o.a-iokcV] 
from the fitting results, we can determine the value of i?in, 
black hole mass, and also T[o.3-iokcV]/-^E, where Le is the 
Eddington luminosity (= 1.5 x 10^®(M/MQ)erg s^^). 
3.2.2. Problem of DBB+PL model fitting 

By using the equations above and the values of Tin and 
/[0.3-iokcV] from the fitting results, we obtained, for the 
case of NGC 5204 X-1, i?in = 2.54 x 10^{cosi)-^^^ km and 
the derived black hole mass of 287/3-1 (cosi)-i/2M0. We 
also obtained T[o.3_iokcV] /Le = 0.05/3(cosi)-i/2^ In the 
same way, we obtained relatively large i?in and hence large 
black hole mass exceeding 100 Mq for other sources (see 
Table 2). 

To summarize, we found through the fitting with the 




Fig. 1. Continued. 



DBB+PL model that the derived black hole masses by far 
exceed 100 Mq and are, hence, within the IMBH regime. 
Thus, fitting the data with the conventional model sup- 
ports the idea of sub-critical accretion onto the IMBHs. 

However, we here point out a serious problem in this 
interpretation. That is, the spectral decomposition shows 
that the power-law component predominates in all energy 
region in question (see Fig. 2 for the case of NGC 5204 X- 
1). In other words, the disk component only gives a minor 
contribution in the entire energy ranges for all the sources. 
Therefore, one is never allowed to use equations (l)-(4), 
which are derived on the assumption that the radiation is 
for 100 % from the DBB component. The black hole mass 
values derived above cannot be so reliable. 

3.3. p-free model 

Next, we tried to fit the same data set by the simplified 
super-critical disk model alone; i.e., the p-free disk model. 
3.3.1. Fitting results 

Firstly, we tried the case of NGC 5204 X-1. Surprisingly 
we also obtained a good fit with this model, although only 
a single component is considered (see Fig.3). The reduced 
chi-squared is 1.12. More importantly, we obtained p = 



0.50 ±0.03, that is just the value predicted by the theory 
of super-critical accretion disk (slim disk; Wang, Zhou 
1999; Watarai, Fukue 1999). We also obtained high inner 
disk temperature, fcTjn = 2.54 ±0.34, as is expected from 
large spectral hardening factors, which the slim disk model 
predicts. 

NGC 5204 X-1 is not an exception (see Table 1). We 
also found a good fit with the p-free model for the cases of 
NGC 4559 X-7 and X-10. The p values of 0.50 ± 0.03 (X- 
7) and 0.51 ±0.06 (X-10) and the inner disk temperatures 
of 1.84 ±0.16 keV (X-7) and 3.16 ±0.71 keV (X-10) are 
again consistent with the slim disk model. 

As for the cases of NGC 1313 X-2, fitting with the p- 
free model only also gave a good fit. We obtained, kT^ = 
2.00 ± 0.34 keV, and the p- value is 0.50 ± 0.05. 
3.3.2. Interpretation of p-free model 

It is clear from the results of the fitting with p- 
free model that the obtained p value is around 0.5, 
indicating super-critical accretion fiow (slim disk) for 
all the data. Applying similar interpretation as in 
the case of fitting with DBB+PL, we obtained, = 
23.7(cosi)~^/^(K/1.7)^ km, the derived black hole mass 
of 2.67p-\cosi)-^/^iK/1.7fMQ and £[o.3-iokcV]/iE = 
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Fig. 2. Spectral decompositions of the fits for NGC 5204 X-1. The thinner lines represent, from bottom to top, disk component, 
power-law component, and total component, while the thicker line represents p-free model 



5.29/3(cosi)-i/2(K/i.7)-2 foj. NGC 5204 X-1. 

Caution is needed here, however, regarding the inner 
edge of the disks in the shm disk regimes and the plau- 
sible spectral hardening factor values (k). Through the 
careful numerical integration, Watarai et al. (2000) found 
that the inner edge of the disk is substantially smaller, 
i?in ~ Rs (i-e., P ~ 1/3), in the slim disk regime, compared 
with that in sub-critical (standard disk) regime, i?in 3i?s 
(/3 ^ 1), even though a central black hole is not rotating. 
This is because enhanced mass accretion flow with high 
density finally fills up an empty zone inside 3i?s , which ex- 
ists in the sub-critical accretion regimes. We then expect 
significant blackbody emission from inside 3i?g . ft is thus 
preferable to set (3^^1/3, instead of /3= 1. However, radia- 
tion from the vicinity of black holes should be attenuated 
because of intense gravitational redshifts and transverse 
Doppler effect, and since the effective temperature cor- 
rected for the rclativistic effect has a peak at around 2i?s , 
/? ~ 2/3 will be a more reasonable estimate (Kawaguchi 
2003). 

As for the spectral hardening factor, Shimura and 
Takahara (1995) and Zampieri et al. (2001) found k « 1.7 
for stellar-mass black holes, and it increases with increas- 
ing M . In the slim disk, the disk temperature gets higher 
and the density becomes smaller as M increases, there- 
fore the electron scattering becomes more dominant and 
the spectral hardening factor becomes much higher than 
- 1.7 (Kawaguchi 2003). Watarai and Mineshige (2003) 
found from the comparison between the limit-cycle insta- 
bility theory (Honma et al. 1991) and the observation of 



GRS1915+105 that the spectral hardening factor close to 
the Eddington luminosity should be about 3. 

The mass estimated from the p-frec models is propor- 
tional to K^, so if we take a larger k, say, 3, the inferred 
mass will be tripled compared to the case of standard 
disks {k ^ 1.7). Consequently, adoping /3 2/3, our mass 
estimate is 12.5(cosi)~^/^(K/3.O)^A'/0 and the Eddington 
ratio is L[o.3-iokcV]/iE l.l(cosi)"^/^(K/3.0)-2 for NGC 
5204 X-1. In the same manner, we find smaller black hole 
masses and larger Eddington ratios for other sources (see 
Table 2). The derived black hole mass is well within the 
stellar mass black hole regime. 



3.4- Fitting Results with Kawaguchi's model 

Finally, we present the results of the data fitting with 
the slim disk model developed by Kawaguchi (2003) (see 
Foschini et al. 2005 for the case of ULX M33 X-8). He 
developed several models, but we chose the model in which 
the effects of electron scattering and relativistic correction 
are included. Calculation was made and fitting was tried 
for several different values of the viscous parameter, a. 
We found a = 0.1 gives most satisfactory fits, and will 
show only the results with a = 0.1 in this paper. Source 
distance was fixed, and the face-on geometry assumed. 
Model parameters are only mass and mass accretion rate 
(in the unit of £Edd/c^). We summarize the fitting results 
using this model in Table 3. 

The best fit for NGC 5204 X-1 was obtained with re- 
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Table 1. Fitting results with DBB+PL and p-frec model. A'^h is an absorption column external to our Galaxy in unit of 10^^ atom 
cm~^. A Galactic absorption column is assumed as in the text. fcTin is the inner disk temperature in keV. The normalization for 
DBB and p-free is defined by [{iJin/km)/(D/10kpc)]'^ cosi, where is the inner disk radius, D the distance to the source, and i 
the angle of the disk. The unit of power-law normalization is photon/s/cm2/keV at 1 keV. /x is the observed flux in units of erg 
cm~-^ (0.3 — 10 keV) that depends on the model. 



Parameter NGC 5204 X-1 



NGC 4559 X-7 



NGC 4559 X-10 



NGC 1313 X-2 



DBB+PL 





0.51±0.11 


1.58±0.19 


0.98 ±0.22 


1.98 ±0.27 


kT 


0.25 ±0.03 


0.17±0.09 


0.48 ±0.25 


0.27 ±0.04 


Norm. 


3.7 ±2.6 


46.9 ±30.3 


0.04±0.11 


4.5 ±4.0 


r 


1.92 ±0.06 


2.17 ±0.05 


1.94 ±0.13 


2.01 ±0.13 


Norm. 


(3.02 ± 0.48) xlO"'' 


(2.10±0.24)xl0-4 


(2.01 ± 0.83) X 10-4 


(2.02 ± 0.73) xlO"'' 


/x 


1.66 X 10^12 


8.36 X 10"" 


9.42 X 10"" 


8.63 X 10"" 


xVdof 


1.05 


0.96 


0.93 


0.98 


p-free 












0.34 ±0.01 


0.60 ±0.06 


0.91 ±0.08 


1.42 ±0.14 




2.54 ±0.34 


1.84 ±0.16 


3.16 ±0.71 


2.03 ±0.35 


P 


0.50 ±0.03 


0.50 ±0.03 


0.51 ±0.06 


0.50 ±0.05 


Norm. 


(5.04±6.11)xl0-4 


(1.04±0.85)xl0-3 


(1.51±3.13)xl0-'' 


(8.13±12.9)xl0-4 


/x 


1.55 X 10-12 


7.96 X 10"" 


9.14 X 10"" 


8.08 X 10-" 


xVdof 


1.12 


1.19 


0.93 


1.15 



duced chi-squared of 1.18, black hole mass of I8M0 and 
mass accretion rate of 27A/Edd, where A/Edd = L^dd/c^ 
^. All the analyzed data of ULXs can be fitted with 
Kawaguchi's model and the derived black hole masses are 
less than 100 Mq (sec Table 3). This is not so surprising 
for the data showing the slim disk signatures, because the 
basic spectral shapes of his model are similar to those of 
the p-free model with p ~ 0.5 (as long as M 3> Le/c^)- 



4. Discussion 

4-1. Accretion Disk Models and Black Hole Mass 

In the present study, we re-examined the XMM-Newton 
data of ULXs through the conventional and more sophis- 
ticated spectral models constructed based on the theoreti- 
cal study of super-critical accretion flow. We would like to 
comment on the results from fitting with DBB+PL and 
p-free model first. It is obvious from the fitting results 
in other papers, e.g. Roberts et al. (2005), Cropper et al. 
(2004), Miller et al. (2003,2004), that the DBB+PL model 
provides a good fit to several ULXs data, and the 'cool' 



Here, 16LEdd/c^ is the critical accretion rate, which would pro- 
duce Eddington luminosity in the case of a classical radiative 
efficiency (1/16). 



disk temperature obtained from the fitting suggests an ex- 
istence of IMBHs inside the ULXs. However, we should 
note that this is only one option, and we tried another 
option; namely, the p-free model. Remarkably, the p-free 
model also gives good fit to the same data. Moreover, the 
obtained p values, which are ~ 0.5, indicating the super- 
critical accretion (slim disk) model, in stead of standard 
disk (DBB) model where, p = 0.75. The derived black 
hole masses and the Eddington ratios also support this 
conclusion. 

Why can the two distinct models give similarly good 
fits to the data? A trick lies in the spectral slope of the 
slim disk model. For a power-law temperature profile, 
T (X r~P, the entire disk spectra, which is the summation 
of the blackbody emission spectra from various radii, give 
power-law spectral slope in the intermediate energy range 
(below hv < fcTin) with a spectral slope of F^, oc 
(see, e.g., chapter 3 of Kato et al. 1998). Then, a photon 
index is F = (2/p) — 2. That is, for a slim disk with p = 1/2, 
we have F = 2. (Note that for a standard disk with p = 3/4 
we get F = 2/3.) This is a striking coincidence that both of 
the PL model with F = 2 and the slim disk model give the 
same power-law spectra in the moderate energy range. To 
see this more explicitly, we show in Fig. 2 the two fitting 
models with fitting parameters obtained for NGC 5204 X- 
1. We see that both DBB+PL model and slim disk model 
produce power-law decline with a power-law index — 1 
(or photon index, F ~ 2) and are identical at 0.5 ^ 5 keV 
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Table 2. The derived black hole mass value from DBB+PL model and p-frce model in Mq and the ratio between luminosity at 0.3 
- 10 kcV and Eddington luminosity. We assumed /3 ~ 1 for DBB+PL model and /3 ~ 2/3 for p-free model (see text). The spectral 
hardening factor, ft, is assumed to be 1.7 for DBB model, and considered to be 3 for slim disk. 



NGC5204 X-1 NGC4559 X-7 NGC4559 X-10 NGC1313 X-2 



DBB+PL 

M(cosi)i/2 287 889 118 136 

L^cosif/^/LE 0.05 0.03 0.3 0.03 



p-free 

M(cosi)i/2 12.5(f)2 34.6(f)2 12.5(f)2 11.2(f)2 

L{cosiY'^lL^ l-l(t)"' 0.9(f)-2 2.8(f)-2 0.4(f)-2 



Table 3. Fitting results with slim disk model of Kawaguchi (2003). Afjj is an absorption column external to our Galaxy in unit of 
10^^ atom cm~^. A Galactic absorption column is assumed as in the text. M is mass in Mq. M is mass accretion rate (= Z/Edd/c^). 
/x is the observed flux in units of 10^^^ erg cm^^ s^^ (0.3 — 10 keV) that depends on the model. 



Fit Par. NGC5204 X-1 NGC4559 X-7 NGC4559 X-10 NGC1313 X-2 



iVn 0.09 ±0.03 0.45 ±0.03 0.69 ±0.07 1.23 ±0.07 

M 18^2 57j:4 24t^ I2t\ 

M 27tl 100^40 141^ 

/x 1.6 0.79 0.91 0.8 

ifo.3-iokcV]/iE 0.81 0.52 1.42 0.36 

xVdof 1.18 1.20 0.93 1.17 



energy range. 

In the present paper, it was hardly possible from model 
fitting to judge if the DBB+power-law model or the slim 
disk modesl is preferred, being limited by photon statis- 
tics. It is, however, in principle possible to distinguish the 
two models by measuring their spectra very accurately. So 
far, XMM-Newton has been able to do this only for M82 
X-1, the brightest ULX (Okajima et al. 2006). Okajima 
et al. (2006) found that the energy spectrum of M82 X-1 
above ^ 3 keV has a curvature just between the power- 
law and DBB model, and is represented pretty well with 
p-free model of 0.61. This is considered to be strong ev- 
idence of slim disk in ULXs. Next generation hard X-ray 
satellites with higher sensitivity and imaging capability 
are expected to precisely measure the energy spectra of 
most ULXs, so that black hole parameters are consrained 
by applying precise slim disk model spectra. 

4.. 2. Luminosity- Temperature diagram 

Then, which is more likely the case, sub-critical accre- 
tion onto intermediate-mass black holes or super-critical 
accretion onto stellar-mass black holes? 

We pointed out a serious problem inherited from the 
DBB+PL model. The spectral fitting with DBB+PL 
model shows that the power-law component predominates 
in all energy regions (Fig. 2). In other words, disk con- 
tribution is very small compared to the power-law contri- 
bution and therefore the black hole masses derived above 



might not be so reliable, because the standard disk rela- 
tions cannot be used here. 

On the other hand, a good fit with p-free model alone 
implies that the disk contributes in all energy region, 
therefore, the black hole masses derived from the fitting 
with p-free model can be reliable. For these reasons, we 
can safely conclude that the latter model is more likely 
the case to explain the physics of ULXs. 

We summarize our results in Fig. 3. This is the 
temperature-luminosity diagram (cf. Makishima et al. 
2000; cf. Miller et al. 2004) with the lines of constant 
black-hole masses. We plot our ULXs fitting results with 
DBB+PL and p-free model together with several best- 
studied stellar-mass black hole candidates (BHCs). We 
can see that the low temperatures obtained from the fit- 
ting using the DBB+PL model are shifted to higher tem- 
peratures when fitted with the p-free model. This implies 
that the black hole mass values in the range of IMBHs 
obtained from the fitting with the DBB+PL model can 
be shifted to the stellar mass range by the p-free model; 
that is, still there is no evidence of IMBHs so far, at least 
in four ULXs investigated in this paper, due to the minor 
disk contribution in DBB+PL model as mentioned before. 

It then naturally follows that all the fitting results which 
gave low disk temperatures and photon indices of F ~ 2 
should be re-examined, since those data are very likely to 
be fitted with the slim disk (p-free) model equally well. 
Obviously, the latter model gives higher disk tempera- 
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Fig. 3. Relation between the luminosity at 0.3 - 10 keV (i[o.3— lOkcV] ) ^^'^ ^in for ULXs and stellar mass black hole candidates 
(BHCs). BHCs references: Ebisawa 1991 (LMC X-1), Treves et al. 1988 (LMC X-3), Mendez et al.l998 (GRO J1655-40), and 
Belloni et al. 1997 (GRS 1915+105). Note that in this figure L]-,oi means the 0.3 - 10 keV luminosity for ULXs and the bolometric 
luminosity for Galactic black hole candidates, both assuming the face-on geometry. The ranges of Tin in the p-free model box 
represent the values via DBB model (lowest end) and p-free model (highest end). 



tures and hence lower black hole masses. It is, in this 
respect, interesting to note the statistical study of ULXs 
by Winter et al. (2006) who examined distribution of 
photon indices and found a concentration at F ~ 2 for the 
high-state ULXs, those ULXs whose spectra can well be 
fitted with the combined DBB+PL model. Some of them 
may be in the slim disk regimes, rather than in the high 
state. 

5. Conclusion 

We have investigated the XMM-Newton EPIC spectra 
of the four ULXs which were claimed to be strong candi- 
dates of IMBHs by several authors. We found that these 
spectra are successfully fitted with slim disk models, and 
the derived masses arc ~11 to 60 Mq depending on dif- 
ferent assumptions. We have not seen the evidence of 
the "intermediate mass" black holes having several hun- 
dreds of Mq . When the mass accretion rate is close to or 
exceeding the critical rates, the slim disk is theoretically 
predicted. We suggest that ULXs are stellar mass black 
holes, at most several tens of Mq under super-critical ac- 
cretion rates shining at super-Eddington luminosities. 
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